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We utilized a simple, robust, first principle method, based on basis set optimization with the 
BZW-EF method, to study the electronic and related properties of transition metal mono-nitrides: 
ScN and YN. We solved the KS system of equations self-consistently within the linear combination 
of atomic orbitals (LCAO) formalism. It is shown that the band gap and low energy conduction 
bands, as well as elastic and structural properties, can be calculated with a reasonable accuracy 
when the LCAO formalism is used to obtain an optimal basis. Our calculated, indirect electronic 
band gap (E^ x ) is 0.79 (LDA) and 0.88 eV (GGA) for ScN. In the case of YN, we predict an 
indirect band gap (Eg~ x ) of 1.09 (LDA) and 1.15 eV (GGA). We also calculated the equilibrium 
lattice constants, the bulk moduli (B D ), effective masses, and elastic constants for both systems. Our 
calculated values are in excellent agreement with experimental ones where the latter are available. 

PACS numbers: 77.84.Bw, 62.20.-x, 71.20.Be, 71.20.-b, 71.15.Mb 



I. INTRODUCTION 

There has been great interest in transition-metal ni- 
trides in the past several decades [H3]. This interest 
is fueled by their many potential, technological applica- 
tions, including high hardness, [5] high temperature sta- 
bility [Sj, mechanical strength (THH], magnetic [TU], and 
electronic properties that vary from semiconducting to 
metallic phases [TlTf2"2"] . These unique properties make 
XN useful refractory materials and hard coatings for cut- 
ting tools [TTHT3] . 

The aforementioned properties of these materials moti- 
vated a variety of experimental [BUS] \TT \ [L2 l H2 | fT5 | I22H27] 
as well as theoretical studies [3 [3 M H3 HI [MM (23- 
135] with the latter utilizing computational techniques of 
varying sophistication, ranging from the tight binding 
and the empirical pseudopotential methods to ab-initio 
density functional theory (DFT). 

Even with the vast number of experimental and com- 
putational studies already available for XN, a satisfactory 
description of their electronic, transport (effective mass), 
elastic, and structural properties is still an area of active 
research. Experimentally, ScN is known to be a semicon- 
ductor with a band gap in the range of 0.90 ± 0.1 eV to 
1.32 ± 0.3 eV [HEI, while most DFT calculations utiliz- 
ing local density approximation (LDA) and generalized 
gradient approximation (GGA) potentials found it to be 
a metal [31 IT51 13111 I31)H39"] . Recent Green's function quasi- 
particle [321 HQ] , exact exchange [36l EI] and screened ex- 
change [T5] 22] calculations have reported that ScN is a 
semiconductor with an indirect gap (E^ _x ) in the range 
0.54 - 1.70 eV. For YN, few theoretical calculations utiliz- 
ing various forms of DFT potentials have reported that it 
is semi-metallic or semiconductor with an indirect band 
gap (E£" x ) of 0.2 -0.3eV HHi], 0.54 eV [35], 
and 0.80 eV [45]. We are not aware of any experiments 
reporting the indirect band gap of YN. 



Theoretical computations have had difficulties in pre- 
dicting the correct band gap energy and other related 
electronic properties of XN from first principle. In- 
deed, the "band gap problem" is decades old. Several 
approaches to solve it have been proposed with signifi- 
cant successes. Density functional theory plus additional 
Couloumb interactions (DFT+U) formalism [4"6TPI9"] has 
had good successes in obtaining correct energy bands and 
gaps of materials, but can only be applied to correlated 
and localized electrons, e.g., 3d or 4/ in transition and 
rare-earth oxides. Hybrid functionals }5T)H52~] have also 
been used in attempts to improve on the calculated en- 
ergy bands and band gaps of materials. This approach 
involves a range separation of the exchange energy into 
some fraction of nonlocal Hartree-Fock exchange poten- 
tial and a fraction of local spin density approximation 
(LSDA) or generalized gradient approximation (GGA) 
exchange potential. We should note that this range sepa- 
ration is not universal. There is always a range separation 
parameter uj which varies between and oo. While it is 
reasonably clear that there exists a value of u that gives 
the correct gap for a given system, this u) is not universal 
as its value is always adjusted from one system to another 
[S3 El]- For example, in HSE06 [53 [55], a; = O.lla^ 1 
(do is the Bohr radius) and in Perdew-Burke-Ernzcrhof 
(PBEh) global hybrid [5S], it is 25 % short-range exact 
exchange and 75 % short-range PBE exchange. Even 
though the HSE functional, in most cases, accurately re- 
produces the optical gap in semiconductors, it severely 
underestimates the gap in insulators [53] [57] and the band 
width in metallic systems is generally too large [531 1571 - 
[5J], The Engel and Vosko [50] (EV) GGA and the Tran 
and Blaha [61] modified Becke- Johnson (TB-mBJ) po- 
tentials have also provided some improvements to the 
calculated band gap of materials. For TB-mBJ, while 
the band gaps are considerably improved, the effective 
masses are severely underestimated. [59 In the case of 
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the EV-GGA potential, the equilibrium lattice constants 
are far too large as compared to experiment and, as such, 
leads to an unsatisfactory total energy j62j [63] . 

The theoretical underestimations of band gaps and 
other energy eigenvalues have been ascribed to the in- 
adequacies of density functional potentials for the de- 
scription of ground state electronic properties of XN . 
Also, other methods [5U [f>5] that entirely go beyond the 
density functional theory (DFT) do not obtain the cor- 
rect band gap values of most semiconductors without ad- 
justment or fitting parameters |29j . This unsatisfactory 
situation is a key motivation for our work. In light of 
the promising technological properties of these materi- 
als, parameter-free computations could aid significantly 
in the design and fabrication of devices. To this end, we 
have also investigated the elastic properties of ScN and 
YN within our parameter-free method. In YN, we pre- 
dict several of these properties for which there are yet no 
experimental data. 

In this paper, we utilized a simple and robust approach 
based on basis set optimization. The rest of this paper 
is organized as follows. After this introduction in Sec- 
tion |lj the computational method and details used in 
our work are given in Section [IT] Section |III| shows our 
computed results. Calculated electronic structures are 



given in Subsection III A The results for the chemical 
bonding and structural properties are presented and dis- 
cussed in Subsections IIII Bl and IIII CI We will conclude 



with a summary in Section IV 



II. 



METHOD AND COMPUTATIONAL 
DETAILS 

A. Method 



We utilized the electronic structure package developed 
at Ames Laboratory of the U.S. Department of Energy 
(DOE), Ames, Iowa [55]. For the LDA computations, we 
used the Ceperley and Alder DFT exchange-correlation 
contribution [67] as parameterized by Vosko-Wilk-Nusair 
[55] , We refer to it as the CA-VWN potential. The 
GGA calculations were carried out using the Ceperley 
and Alder DFT exchange correlation contribution [67] as 
parameterized by Perdew and Wang [B5]. We refer to it 
as the CA-PW potential. 

We employed the linear combination of atomic orbitals 
(LCAO) approach where an unknown wave function, for 
the solid state, is expressed as a linear combination of 
atomic orbitals of the individual atoms in the system. 
The radial parts of these orbitals are generally exponen- 
tial or Gaussian functions resulting from self consistent 
calculations of energy levels of the atomic or ionic species 
that are present in the solid under study. We use Gaus- 
sian functions and refer to that rendition of LCAO as the 
linear combination of Gaussian orbitals (LCGO). 

In addition to the use of DFT potentials and of the 
LCAO formalism, our computational approach rests on 



the Bagayoko, Zhao, and Williams [TDJ, [71] method as 
recently enhanced by the work of Ekuma 72 75] and 
Franklin [76] (BZW-EF). The method heeds the state- 
ment by Kohn and Sham [77] that the equations defining 
LDA have to be solved self-consistently. Upon the selec- 
tion of an exchange correlation potential, these equations 
reduce to (1) the one giving the ground state charge den- 
sity in terms of the wave functions of the occupied states 
and (2) the Kohn-Sham equation. We describe below the 
aforementioned method. 

In the BZW-EF method, we begin our calculations 
with the minimum basis set, the one that is just large 
enough to account for all the electrons in the system un- 
der study (in this case XN). Our self consistent calcula- 
tion with this basis set is followed by another where the 
basis set is augmented with one additional orbital from 
one of the atoms or ions in the system. The compari- 
son of the occupied energies from Calculations I and II 
generally shows that they are different, with those from 
Calculation II generally lower than their corresponding 
ones from I. A third, self consistent calculation is per- 
formed with a basis set that includes that for Calcula- 
tion II plus another orbital from one of the atoms in the 
system. This process of augmenting the basis set and of 
carrying out self consistent calculations continues until 
a calculation, say N, is found to have the same occupied 
energies as calculation (N + 1) that follows it, within 
computational uncertainty of about 50 meV. This con- 
vergence of the occupied energies identifies the basis set 
of Calculation N as the optimal one. The optimal ba- 
sis set is thus the smallest one with which the occupied 
energies verifiably reach their respective minima [731 178j . 
The resulting eigenvalues and the corresponding eigen- 
functions provide the physics description of the system. 

The BZW-EF method has been shown to lead to ac- 
curate ground state properties of many semiconductors: 
c-InN [79], w-InN [8Q], w-CdS [73], c-CdS [78], rutile- 
Ti0 2 [75], SrTi0 3 [51], AlAs [El], GaN, Si, C, Ru0 2 
[7T] , BaTi0 3 [70], and carbon nanotubes [53] , 

In calculations with basis set larger than the optimal 
one, the ground state charge density does not change, nor 
do the Hamiltonian and the eigenvalues of the occupied 
states. Consequently, these calculations do not lower any 
occupied energies (as compared to the results obtained 
with the optimal basis set), even though they generally 
lead to some lower, unoccupied energies [751 [73] by virtue 
of the Rayleigh theorem [5H 155] . The Rayleigh theorem 
can be stated as follows: if one solves an eigenvalue equa- 
tion with two basis sets I and II, with set II larger than I 
and where I is entirely contained in II, then the eigenval- 
ues obtained with set II are lower than or equal to their 
corresponding ones obtained with basis set I (i.e., E^ +1 
< Ef, V i ^N). This theorem explains the reasons that 
some unoccupied energies are lowered when the Kohn- 
Sham equation is solved with basis sets larger than the 
optimal one. Such a lowering of unoccupied energies with 
basis sets larger than the optimal one is fundamentally 
different from the one that occurs before the occupied en- 
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ergies reach their minima. The latter lowering is ascribed 
to physical interactions given that both the charge den- 
sity and the Hamiltonian change, from one calculation to 
the next, before one reaches the optimal basis set, while 
the former one is unphysical. 

Important points about this method include the fact 
that it provides the needed variational freedom for the re- 
structuring of the electronic cloud in the material under 
study, as opposed to the case in isolated atoms or ions. 
Indeed, the methodical increase of the basis set remedies 
possible radial, angular symmetry, and size deficiencies 
by the time the optimal basis set is reached. By changing 
the basis set, i.e., the external input to the charge density 
equation, the method solves self-consistently both this 
equation and the Kohn Sham equation in a way that is 
quintessentially different from input changes limited to 
new expansion coefficients in the iterative process |86j . 
The latter input changes cannot remedy any significant 
deficiency of the basis set, i.e., in terms of radial orbital, 
angular symmetry, or the total number of functions (the 
dimension of the Hamiltonian) . In light of the Rayleigh 
theorem, a deficiency could be over-completeness as well, 
in single trial basis set calculations. The critical impor- 
tance of the BZW-EF method resides in part in the fact 
that, by its very variational derivation, the validity and 
the related physical content of the eigenvalues of the Kohn 
Sham equation directly hinge on the final, self- consistent 
charge density being as close as possible to that of the real 
material under study. 

The thesis work of Williams [71] showed that the 
method may not be needed for the description of materi- 
als established, experimentally, to be metals. At least one 
band crosses the Fermi level in a metal and the lowest- 
laying conduction bands reach their minima at the same 
time as the valence bands do. Consequently, the basis 
set and variational effect [701 EI] that explains the lower- 
ing of unoccupied energies, while the Hamiltonian does 
not change, does not affect these lowest-laying conduc- 
tion bands in metals, even though it lowers some other 
conduction bands. This fact partly explains the early 
success of DFT for the description of metals as opposed 
to semiconductors and insulators. In the initial BZW 
method, we increased the basis set by adding orbitals in 
the order of increasing energies in the atomic or ionic 
species. Recent works by Ekuma 72 7o] and Franklin 
[75] led us to the realization that, for the valence states, 
polarization (i.e., p, d, and / orbitals) have primacy over 
spherical symmetry (s orbital) in diatomic molecules and 
more complex systems with three or more atomic or ionic 
sites. Hence, while still useful, adding orbitals in the or- 
der of increasing energies in atomic or ionic species can 
be relaxed. 



B. Computational Details 

In the ground state, both ScN and YN have the rock- 
salt crystal structure. We utilized room temperature 



experimental lattice constant of 4.501 A [7J I5TH59"] and 
4.837 A [SI] for ScN and YN, respectively. Each Sc (or 
Y) atom is surrounded by six N atoms, thereby providing 
the octahedral environment at the Sc (or Y) site, which 
leads to the splitting of the degenerate d orbital into t2 g 
and e g states. Every Sc (or Y) atom has 12 nearest neigh- 
bors Sc (or Y) atoms and 6 next nearest neighbors Sc (or 
Y) atoms. Preliminary calculations indicated that scan- 
dium and yttrium are closer to Sc 3+ and Y 3+ than to 
the neutral Sc and Y, respectively. Similarly, nitrogen 
is N 3 ~ as opposed to the neutral N. Therefore, we first 
carried out self consistent calculations of the electronic 
properties of X 3+ (X = Sc, Y) and N 3 ~. Atomic or- 
bitals utilized in these calculations, for the valence states, 
are given between parentheses: Sc 3+ (Zs'ipid^sApAd) and 
N 3 - (2s2p3s3p) for ScN and Y 3 + (AsApAdbsbp) and N 3 " 
(2s2p3s3p) for YN. Other atomic states with higher bind- 
ing energies were treated as deep core states. In the op- 
timal basis set for the valence states of ScN, (ApAd) and 
(3p) are unoccupied for Sc 3+ and N 3_ , respectively. For 
YN, the unoccupied orbitals in the optimal basis set for 
the valence states are (bp) and (3p) for Y 3+ and N 3_ , 
respectively. Nevertheless, these orbitals are included in 
the self-consistent LCAO calculations to allow for a re- 
organization of electronic cloud in the solid environment, 
including polarization. 

In the self-consistency calculation, both the potential 
and charge density are also expanded in terms of even 
tempered Gaussian orbitals (Sc: 15, 15, and 13; Y: 15, 
15, and 13; N: 17, 17, and 15 for s, p, and d orbitals, 
respectively). The exponents, a, of the Gaussian basis 
sets range from 0.19 to 10 5 . The charge fitting error 
using the Gaussian functions in the atomic calculation 
was about 10~ 5 . Since the deep core states are fully 
occupied and are inactive chemically in the materials, 
the charge densities of the deep core states are kept the 
same as in the free atoms. However, the core states of 
low binding energy are still allowed to fully relax, along 
with the valence states, in the self consistent calculations. 
The computational error for the valence charge was 1.2 
x 10~ 5 and 5.9 x 10~ 6 per valence electron for ScN and 
YN, respectively. The self consistent potential converged 
to a difference of 10 -5 after about 60 iterations. 

The Brillouin zone (BZ) integration for the charge den- 
sity in the self consistent iterations was based on 28 spe- 
cial k points in the irreducible BZ (IBZ). The energy 
eigenvalues and eigenfunctions were then obtained at 161 
special k points in the IBZ for the band structure. A total 
of 152 weighted k points, chosen along the high symmetry 
lines in the IBZ of ScN and YN, respectively, were used 
to solve for the energy eigenvalues from which the elec- 
tron densities of states (DOS) were calculated using the 
linear, analytical tetrahedron method [90]. The partial 
density of states (pDOS) and the effective charge of each 
atom were calculated using the Mulliken charge analysis 
procedure [91] . We also calculated the equilibrium lattice 
constant (a Q ), the bulk modulus (B ) and the electron ef- 
fective masses for different directions with respect to the 
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r point. For the calculation of the equilibrium lattice 
constant, we utilized the Murnaghan equation of state 
[92), [93] . By applying an appropriate set of strains to 
the undeformed unit cell lattice, we calculated the elas- 
tic constants from the resulting change in total energy 
on the deformation using the strain-energy method. For 
a typical cubic crystal, three independent elastic moduli 
are of importance; they are usually denoted as Cn, C12, 
and C44. 



III. RESULTS AND DISCUSSION 

The results of the electronic structure computations 
are given in Figs.[T]to [2] Figure [3] shows the contour plot 
of the distribution of the electron charge density while 
Figure [4] depicts the total energy curves of XN. We dis- 
cuss the electronic structure and the effective masses in 
|III A| Densities of States, chemical bonding, and electron 
distribution are described in Subsection IHI Bl The struc- 
tural and elastic properties are presented in Subsections 
IITTC1 



The Electronic Structure, Band gap and 
Effective mass 



Figures [T]to [2] exhibit the energy bands, and the re- 
lated total (DOS) and partial (pDOS) densities of states 
of XN. All energies are referred relative to zero energy at 
the top of the valence band (VB). The electronic struc- 
tures of the valence bands, the low energy conduction 
bands, and the band gap determine the most important 
properties of these materials in device applications. 

Our ab-initio method shows that the fundamental gaps 
of both ScN and YN are indirect ones, with the maxi- 
mum of the valence band (VB rnax ) occurring at T and 
the minimum of the conduction band (CB m i n ) at the 
X point (see Figs. [j|a) and [^a) for ScN and YN, re- 
spectively). Our calculated, indirect gaps for ScN, using 
the equilibrium lattice constant, is 0.79 eV and 0.88 eV 
for CA-VNW (LDA) and CA-PW (GGA) potentials, re- 
spectively. The difference in the band gaps and energy 
eigenvalues between CA-VWN and CA-PW may be at- 
tributed to the enhancement factor (s = ) [9~4l |9"5] 
in the GGA functional. The band gap of ScN has been 
reported to suffer from the Burnstein-Moss shift \§T\ 
due to large background carrier concentration. As a con- 
sequence of this effect, several experimental indirect band 
gap values ranging from 0.90 ± 0.1 eV [4 to 1.32 ± 0.3 eV 
[98] have been reported. Our computed band gap values 
are in good agreement with the lower, experimental val- 
ues that result from relatively low carrier concentrations. 
For YN, our computations predict an indirect band gap 
value of 1.09 eV and 1.15 eV for LDA and GGA, respec- 
tively. To our knowledge, no measurement of the YN 
indirect gap has been reported. We expect YN to suffer 



from the same Burnstein-Moss effect as ScN, due to the 
close similarity of their structures. 

Our calculated valence bands for ScN (see. Fig. [ija)) 
and YN (see Fig. [2^ a)) resemble those from other calcu- 
lations available in literature. A major difference is the 
size of the band gap. For both ScN and YN, there is sig- 
nificant N 2p - Xid (X = Sc, Y; and i = 3, 4 for Sc and Y, 
respectively) hybridization in the valence bands. The low 
energy conduction bands up to about 5.24 eV are mainly 
from X ld (see Figs. Qb) and ||b) for ScN and YN, re- 
spectively). The following can further be confirmed from 
the partial density of states (see Figs, [ijc) and [2^c) for 
ScN and YN, respectively). The upper valence bands 
originated from the bonding between N 2p states and an 
admixture of 3d states (for Sc) and Ad states (for Y). 
The lowest conduction bands are mainly the antibonding 
t 2g states with 3d (for Sc) and id (for Y) character, re- 
spectively. The two nonbonding e g bands are higher in 
energy. These bonding characters are in perfect agree- 
ment with the bonding analysis of Harrison and Straub 
[9"9"] . The authors showed that rock-salt XN have 3 p-like 
bonding, 3 d-like antibonding t2 g , and two c?-like non- 
bonding e g bands formed by the hybridization of three 
valence p states of N with the five d states of X. The 
assignment of the bondings is consistent with the partial 
density of states analysis of the screened exchange LDA 
FLAPW calculations of Stampfl and co-workers [T51 [3T] 
and of the exact exchange based quasiparticle calcula- 
tions of Qteish et al. respectively. 

The density of states (DOS) of ScN is shown in Fig. 
[ljb). In the valence bands, we found sharp peaks at 
-2.46 ± 0.1 eV and -2.96 ± 0.1 eV, respectively. A shal- 
low minimum is observed at -3.05 ±0.1 eV, followed by 
a shoulder at -3.2 ± 0.1 eV. A broad peak is located at 
-4.05 ± 0.2 eV. These peaks are formed by a strong hy- 
bridization between N p states and Sc d states with little 
contribution from Sc p and s states, respectively. In the 
conduction bands, the first observed, significant feature 
is a shoulder at about 3.80 ±0.1 eV. It is due to the hy- 
bridization between N p and Sc d states. A sharp peak is 
observed at 4.7 ± 0.1 eV; it is formed from all the states 
of the atomic species except N s states. The density of 
states of YN (see Fig. |2jb)) is similar to that of ScN, ex- 
cept that it has a pronounced shoulder at 0.92 eV below 
the Fermi energy. In the DOS of the valence bands, we 
predicted a sharp peak at -1.72 ±0.1 eV, followed by 
two smaller ones at -2.1 ± 0.2 eV and -2.82 ±0.1 eV, 
respectively. These peaks are mainly of N p states and Y 
d states. In the conduction bands, we observed a small 
peak at 3.6 ± 0.1 eV, which is mainly of N p states and 
Y d states. A sharp peak is observed at 5.94 ± 0.1 eV. 
This peak is of N s, N p, and Y d states character. 

To determine some specific properties of XN relevant to 
transport, we calculated the effective mass in different di- 
rections with respect to the V point. The effective masses 
are calculated from curves that fit the bands in the imme- 
diate vicinity of minima (for electrons) or a maxima (for 
holes). The calculated electron effective masses (in units 
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FIG. 1: (Color online), (a) The calculated band structure of c-ScN as obtained with the optimal basis set (for the equilibrium 
lattice constant of 4.501 A). The Fermi energy (Ef) has been set equal to zero, (b) The calculated density of states (DOS) of 
c-ScN, as obtained from the bands shown in Fig. [TJa). The Fermi energy (Ef) has been set equal to zero, (c) The calculated 
partial density of states (pDOS) of c-ScN, as obtained from the bands shown in Fig.flja). The Fermi energy (Ef) has been set 
equal to zero. 




FIG. 2: (Color online), (a) The calculated band structure of c-YN as obtained with the optimal basis set, at the equilibrium 
lattice constant of 4.827 A. The Fermi energy (Ef) has been set equal to zero, (b) The calculated density of states (DOS) of 
c-YN as obtained from the bands shown in Fig. ^ a) The Fermi energy (Ef) has been set equal to zero, (c) The calculated 
partial density of states (pDOS) of c-YN, as obtained from the bands (one with solid line) shown in Fig. |2ja). The Fermi 
energy (Ef) has been set equal to zero. 



of m ) at the bottom of the conduction band, along the 
T-L, T-X, and T-K directions, respectively, are 0.91-1.66, 
0.78-0.99, and 0.87-1.07 for ScN, while those for YN are 
0.89-1.02, 0.57-0.61, and 0.64-0.68 in the same direction 
as for ScN. 



B. Chemical Bonding and Electron Charge 
Distribution 



Figure [3] shows the contour plots of the distribution of 
the electron charge densities of XN along the (100) plane, 
cutting through the atoms. Away from the atomic cen- 
ters (i.e., in the interstitial region) the electron charge 
density distributions are not spherically symmetric. The 
chemical bonding in ScN and YN appears to be of inter- 
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FIG. 3: (Color online) (a) The contour plot of the electron 
charge density of ScN as obtained with the optimal basis set. 
(b) The contour plot of the electron charge density of YN as 
obtained with the optimal basis set. A n(r) is the variation 
of the electron charge density as a function of distance away 
from an atomic site. A logarithmic scale has been used. 



mediate character between ionic bonding in the noncon- 
ducting calcium nitride or strontium nitride and that of 
the metallic bonding prevailing in the conducting group 
IV transition nitrides. As a consequence of this com- 
plex bonding, there seems to be a cooperative/competing 
bonding mechanisms in XN. As can be seen from Fig. [3j 
the bonding between X and N is covalent, with a signifi- 
cant ionic character. The bond length of Sc - N is 2.25 A 
while that of Y - N is 2.42 A. The experimental Sc - N 
bond length is 2.24 A pJIHE] while the experimental bond 
length, Y - N, is 2.44 A [100] . The covalent bonding in 
ScN can be seen to be stronger than that in YN. This 
difference is understandable in light of the larger Y - N 
bond. 



C. Structural and Elastic Properties 

Total energy versus lattice constant data are as shown 
in Fig. [4] for ScN and YN. The computation was done 
with both the CA-VWN and CA-PW potentials. How- 
ever, the reported data are for CA-PW, as we noted that 
the total energy difference between results from the two 
potentials were small (± 0.12 eV). The data are fitted 
to the Murnaghan equation of state (EOS) for ScN and 
YN, respectively. The fit is very good, indicating that a 
different choice of EOS and/or lattice constant will have 
insignificant effect on the results. 

For ScN, the calculated equilibrium lattice constant is 
4.501 A, the bulk modulus (B a ) is 200.19 - 216.73 GPa, 
and its pressure derivative (B' ) is 3.32 GPa. Experi- 
mental values are 4.501 A S [55], 182 ± 40 GPa [7], 
and 3.31 - 3.89 GPa for the lattice constant, bulk mod- 
ulus, and pressure derivative, respectively. For YN, our 
computations show that the equilibrium lattice constant 
is 4.827 A, the bulk modulus (B ) is 166.99 GPa. Our 
calculated equilibrium lattice constant of 4.827 A is prac- 
tically the same as the experimental value of 4.837 A; it 
is just 0.21 % smaller than the experimental value. 
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FIG. 4: (Color online) The total energy (Ex) per unit cell as 
a function of lattice constants (A) for c-XN. Our calculated 
equilibrium lattice constant for ScN is 4.501 A (exactly the 
same as the experimental value). The calculated equilibrium 
lattice constant of YN is 4.827 A, 0.21 % smaller than the 
experimental value of 4.837 A. 



Using the elastic constant relations as in Ref . [101] , we 
calculated the independent elastic constants Cn, C12, 
and C44. The calculated independent elastic constants is 
used to determine the compliance constants Sy (Sn, S12, 
and S44) [101] 1102] . The Poisson ratios v are measure of 
materials tendency to react to applied strain. For a cubic 
system (averaged over transverse directions in the [100], 
[110], and [111] directions), they are given by: [1031 1104] 

^<ioo) = - sE, "(no) = " Sll+ 2s 12+ s 44 . and ^m) = - 

S 2s^ 3 +2s7 2 +sii ■ The isotropic shear modulus G [US] can 
be estimated from the average of the shear modulus in 
the Voigt approximation (Gy — l/5(Cn - Ci 2 + 3C44)) 
[106] and Reuss approximation (Gr = 4c^4+3^Ci7-Cil) ) • 
[T07] It is given as G = (G v + G R )/2. Note that the 
Voigt and Reus approximations represent the upper and 
lower bounds of the isotropic shear modulus. In a typi- 
cal cubic system, the isotropic bulk modulus K, can be 
calculated using the relation K — \{G\\ + 2C\?). Using 
these quantites, we obtain the isotropic Young's modulus, 



E 



9KG 
3K+G 



and the isotropic Poisson ratio, 



3K-2G 



2(3K+G) 

~ 5(^(100) +^(110) +^<iii))- The equilibrium elastic con- 
stants, compliance constants (Sn, S12, and S44), the elas- 
tic moduli corresponding to the determined stiffness con- 
stants and other elastic parameters are shown in Table [T] 
The calculated, independent elastic constants for ScN 
are Cn = 452.55 GPa, C 12 = 98.82 GPa, and C 44 = 
184.60 GPa. Using the Zener's anistropy index (ra- 
tio) [ISElEl], V = 2C 44 /(Cii - C12), we quantify the 



7 



TABLE I: Calculated equilibrium lattice parameters a , stiffness, compliance and moduli constants, Young's modulus, Poisson 
ratio, shear modulus, and isotropic bulk modulus of ScN and YN. 

Nitride a D (A) C tJ (GPa) v m) Cy (10~ 3 GPa' 1 ) E fafc) (GPa) G (GPa) K 

Material Cn C12 C44 ^(100) ^(110) ^(111) v_ S11 S12 S44 E( 10 o) E( 11(1 ) E( ln ) E Gy Gr G 

ScN 4.50 453 99 185 0.18 0.17 0.17 0.17 2.40 -0.43 5.42 417 428 431 426 182 181 181 217 

YN 4.83 405 93 145 0.19 0.20 0.20 0.20 2.70 -0.50 6.88 371 355 350 358 150 149 150 197 



anistropy as 1.04. This quantity contains the same infor- 
mation as the ratio of the directional Young's modulus 
(E^, £ is direction) E^ m )/E( 100 ) , suggesting that ScN is 
stiffer in the (111) than (100) direction. The elastic mod- 
uli ratio E (1 oo>:E( U o>:E(iii) is 1.00:1.03:1.03. For YN, 
the calculated, independent elastic constants are: Cn = 
405 GPa, C12 = 92.70 GPa, and C 44 = 145.28 GPa. The 
anisotropy ratio i] is 0.93. The elastic moduli ratio is 
E ( ioo>:E (UQ> :E (m) = 1.0:0.96:0.94. 

For ScN, E = 425.68 GPa, G = 181.47 GPa, and v 
= 0.17. Our computed results are comparable to experi- 
mental results of Gall et al. [7j, Moram et al. |110j . and 
Gall et al. [89 on epitaxial films of ScN. These authors 
reported values of E in the range 270 ± 25 to 388 ± 
20 GPa and v value of 0.15, 0.188 ± 0.002 and 0.20 ± 
0.04. Our calculated and measured v values of ScN are 
in the same range. Our calculated value of E is close 
to the upper bound, given above, of the measured val- 
ues for the films. The slight differences between some 
of our computed elastic constants and experimental ones 
may be due to the fact that these experiments were done 
on films of ScN. We should note that growth conditions 
and film thickness [7] are known to affect significantly the 
elastic constants, as is evident in the works of Gall et al. 
(7j, Gall et al. [5S], and Moram et al. (110] which showed 
a wide range of values. For YN, E = 357.94 GPa, G = 
149.53 GPa, and v = 0.20. We note that our computed 
Poisson ratios, v are comparable to known Poisson ratios 
for other transition-metal nitrides that are in the range 
0.19 - 0.22 (Ml EU E2]. We are not aware of any ex- 
perimental measurements of E, G, and v for YN; hence, 
our calculated results are predictions. In both ScN and 
YN, our computed isotropic bulk modulus, K is in basic 
agreement with the bulk modulus we obtained for the 
equilibrium structural optimization. 

Our calculated elastic stiffness constants, obey the sta- 
bility conditions for a cubic crystal system, i.e., B Q = 



(Cn + 2Ci 2 )/3 > 0, G' = (Cn - C 12 )/2 > 0, and C44 > 
[1 13] . The calculated 77 values are both close to 1, im- 
plying that the elastic behavior is almost isotropic. The 
Cauchy pressure (C12 - C44), in both ScN and YN, is 
significantly less than 0. It corresponds to a directional 
bonding [TH I114| 1115] , leading to large charge transfers 
from cations to anions as has been observed experimen- 
tally in TiN [TT5] . 

IV. SUMMARY AND CONCLUSION 

We have performed first principle calculations of elec- 
tronic and related properties of ScN and YN, based on 
basis set optimization with the BZW-EF method. Accu- 
rate, band gaps, structural parameters, lattice constants, 
bulk moduli, and elastic constants were calculated. We 
have provided detailed analyses of the electronic and re- 
lated properties of both ScN and YN. Our computed 
properties are in good agreement with experiment. In 
particular, we predict various electronic (indirect band 
gap value), effective masses, elastic and structural prop- 
erties of YN. It is our hope that the present results will 
motivate further experimental and theoretical studies of 
the properties of these materials, with emphasis on mea- 
surements of electronic and elastic properties of YN. 
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